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Spin excitation spectrum of the double-exchange model is studied in the presence of random- 
ness. Spin wave approximation in the ground state shows that the randomness significantly 
modifies the spectrum from the cosine-like one in the pure system to that with anomalies such 
as broadening, anti-crossing and gap opening. The origin of anomalies is speculated to be modu- 
lation of effective ferromagnetic coupling by the Friedel oscillation. These anomalies qualitatively 
reproduce the spin excitation spectrum in colossal magnetoresistance manganites whose Curie 
temperatures are relatively low. Our results suggest that randomness control is an important 
notion to understand effects of the A-site substitution which has previously been understood as 
the bandwidth control. 
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One of the most important issues of colossal magne- 
toresistance (CMR) manganites AMnOs is to understand 
effects of the A-site substitution.EP Generally, the A-site 
substitution has two different aspects. One is the car- 
rier doping which is the substitution by A ions with 
different valences. The other is the ionic radius modi- 
fication which is the A-site substitution with ions having 
the same valence but different ionic radius, fn the lat- 
ter, the averaged radius of A ions (rA) affects length and 
angle of Mn-O-Mn bonds. This leads to a change of ef- 
fective transfer energies between Mn sites. Therefore, 
it is widely accepted that the ionic radius modification 
gives a bandwidth control. Owing to a wide potential of 
these A-site substitutions, manganese oxides exhibit rich 
physics including the CMR phenomena. 

Among various compounds AMnOs, recent develop- 
ments in both experiments and theories have been re- 
vealed that Lai_a;Sra;Mn03 (LSMO) at f/3 is a 
canonical double-exchange (DE) systems.B© The DE 
modeleP can explain magnetic and transport properties 
in this compound quantitatively. For instance, the exper- 
imental values of Curie temperature Tq agree well with 
the theoretical estimate.cP Thus, the LSMO compound 
provides a well-established starting point to examine ef- 
fects of the A-site substitutions. 

Starting from LSMO as a reference, we focus on the 
ionic radius control by the A-site substitutions. This 
causes a systematic change in transport and magnetic 
properties. For instance, the value of Tq decreases sys- 
tematically for compounds with smaller (rA) and hence 
with 'narrower-bandwidth' such as Pri_j,Srj,Mn03 and 
Lai_a;Caa;Mn03 (LCMO). Hereafter, we call these A-site 
substituted compounds 'low-Tc manganites' whereas we 
call the compounds such as LSMO, which show relatively 
high Tc, 'high-Tc manganites'. An empirical relation be- 
tween Tc and (rA) has been proposed experimentally.cF 



The systematic change has been considered to be pri- 
marily due to the bandwidth change by the A-site sub- 
stitution. 

However, the picture of the bandwidth control has 
been doubted by some experimental results. One is a 
quantitative comparison between Tc and the bandwidth 
estimated from (rA). The change of Tc is much larger 
than that of the bandwidth. For instance, from LSMO 
to LCMO, Tc decreases about 30% while the estimated 
bandwidth decreases less than 2%iS> This large decrease 
of Tc cannot be explained by the bandwidth control in 
the DE theory in which Tc should scale to the band- 
width. An extra degree of freedom as a hidden parameter 
seems to be necessary to explain such a change. 

fn this Letter, we examine a role of the randomness 
in the A-site substitution, fn general, oxides are known 
to be far from perfect crystals. Moreover, in mangan- 
ites, randomness is inevitably introduced by the A-site 
substitution since the compounds are solid solutions of 
different A ions. Random distribution of A ions induces 
a charge and structural disorder in the system, which 
introduces randomness in the electronic Hamiltonian. 

fmportance of the randomness effects has been dis- 
cussed both experimentally and theoretically. One is 
that Tc is much affected not only by the averaged ra- 
dius of A ions but also by the standard deviation of the 
ionic radii. A linear decrease of Tc asLa function of the 
standard deviation has been reported.LF Another exper- 
imental fact is a correlation between the residual resis- 
tivity and Tc. The A-site substitution from high-Tc to 
low-Tc manganites causes a systematic incmase of the 
residual resistivity even in single crystals BlP These re- 
sults clearly indicate that it is necessary to take account 
of effects of the randomness in the A-site substitution. 
Recently, effects q£ tj3£_cajidomness have been discussed 
in the DE modeI.Ej'EiESEj' A significant decrease of Tc 
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by the randomness has been predictedJi£l'til'lliHl§lla) 

Let us now consider another important result of the 
ionic radius modification which is the systematic change 
in the spin excitation spectrum. High^Tc manganites 
show a cosine-hke spin-wave dispersiorO-* which is well 
described by the DE theory quantitatively.Ea^ On the 
other hand, in low-T^ manganites, the spin excitation de- 
viates from this behavior considerably, and shows some 
anomalies such as-aafiiiiiing, broadening, anti-crossing 
and gap-opening.E3E3E^'E3) The DE theory also fails to 
explain these anomalies in the spin dynamics. 

In this Letter, we study effects of the randomness 
on the spin dynamics in the DE model. The spin ex- 
citation is calculated in the ground state by the spin 
wave approximation. We find clear deviation from the 
cosine-like spin-wave dispersion in the pure DE model 
and some anomalies in the spin excitation spectrum. 
Through comparisons with the experimental results in 
low-Tc manganites, our results strongly suggest that the 
randomness gives an important effect of the A-site sub- 
stitution. 

As a model for the A-site substituted manganites, we 
consider the DE model with randomness. The Hamilto- 
nian is written by 

Jh V- 



H 



h.c.)--^ 
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where the first two terms denote the DE HamiltonianI 
which consists of the nearest-neighbor hopping and the 
Hund's-rule coupling between itinerant electrons and lo- 
calized spins with a magnitude S. These two terms quan- 
titatively describe physical propcrtixjSi^ ia. high-Tc man- 
ganites including the spin excitation.aEl'ES' The last term 
denotes the on-site randomness. In this work, for sim- 
plicity, we incorporate the randomness by the A-site sub- 
stitution as the on-site random potential although it may 
cause various other effects, for instance, a randomness in 
transfer integrals. Here we consider the binary-alloy type 
randomness, that is, takes ±VFinip/2 in equal proba- 
bility in each site. We confirm that the results are in- 
dependent of the type of randomness qualitatively. Al- 
though it is difficult to determine the actual magnitude 
of the randomness in manganites, the value of Wimp is 
roughly estimated to be the same orderj:if magnitude as 
the bandwidth in low-Tc compounds. EIlS' As an energy 
unit, we use the half bandwidth at the ground state for 
Wimp = = 0, i.e., W 6t = 1 in three dimensions. 

We study the spin excitation of the model (|l|) using the 
spin wave approximation. We consider a perfectly spin- 
polarized ferromagnetic ground state Si ~ (0, 0, S), and 
calculate the one-magnon excitation spectrum within the 
lowest order of 1/S expansion. In the absence of impu- 



rities, the formulation is given in Ref. 18 



Here we extend the method to the case with random- 
ness. For a given configuration of the quenched random- 
ness {e} on a finite size cluster, we explicitly diagonalize 
the Hamiltonian matrix at the ground state. We have 

Hipnaii) = {En - O" Jh )<^«<t ( «) , (2) 

where (pna- is the n-th wavefunction for the spin a elec- 



tron. Since wavefunctions are given in real space, it is 
straightforward to formulate Green's functions of elec- 
trons and spin wave using the real-space representation. 
Electron Green's function is given by 

^ - {En - o- - m) + iv sgnw 



Gij.a{^) 



(3) 



Following Ref. spin wave self-energy is calculated 
from electron's spin polarization function depicted in 
Fig. 1. In the limit of Jn/t oo, which is realistic 
in manganites, we obtain 
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/«T'^nT(*)^"T(jT<^mi(i)<^mi(«)* 
X {En^ - En " Lo) , (4) 



where N is the number of lattice sites and fn] is the 
fermi distribution function for up-spin states. Spin wave 
Green's function is given by Dij{uj) = {ijj—YVij{uj)-\-'nf)^^. 
Spin wave excitations are obtained from the poles of the 
Green's function. Since Hy oc 1/S, positions of poles 
are at uj ^ 0(1/5'). Then, as discussed in Ref. |l8|, spin 
wave excitations within the lowest order of 1/S expan- 
sion can be determined from the static part of the self- 
energy %j{uj^Q). 




Fig. 1. Spin wave self-energy in the lowest order of the 1/5 ex- 
pansion. See Ref. for details. 

The spectral function for spin wave is defined by 
i(q,w) = VijmAj(t^)exp[iq-(r, -rj)]. (5) 

Using the eigenvalues Coi and eigenvectors ipi{j) of H^, 

^n,,(c^ = 0)-V3/(j) (6) 

j 

the matrix inversion of {uj — Hy ) is obtained. We have 
i(q,c.)-^i,(q)5(c.-ci0, (7) 



^'(9) = ■j;^\^Mi)'^'^'p{'^qrj) 



(8) 



Finally, we average A{q, lo) for random configurations 
of the potential {e} to obtain the spectral function 
A{q,uj). 

Figure 2 shows the spectral function for spin excita- 
tions at a; = 0.3 in three dimensions. Here, x is the hole 
density in model (|l|). System size A^isl6xl6xl6 with 
the periodic boundary conditions and the random aver- 
age is taken for 24 different configurations of the random 
potential for each value of Wmp- For small values of 
Wimp, the spectrum is similar to the cosine-like disper- 
sion for Wimp = (gray curve) although it shows some 
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broadening. However, when the value of Wimp becomes 
moderate, some anomahes appear conspicuously in the 
spectrum. Apart from the large broadening, a remark- 
able feature is a branching which is clearly seen, for in- 
stance, along the F-M and F-R lines for Wimp/W^ — 0.5. 
Near the zone boundary, extra branches with consider- 
able weight appear well below the cosine-like branch. If 
one follows only the extra lower branch, the spin excita- 
tion spectrum appears to show a significant softening. 



pend on impurity strength. From these observations, it 
is clear that the origin of the anomalies is the fermionic 
responses of the itinerant electrons to impurities. 
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Fig. 2. Spin excitation spectra at a; = 0.3 in three dimensions 
for (a) Wi^p/W = 0.2 and (b) W^^p/W = 0.5. The intensity 
is shown by gray scale along the vertical axis. The discreteness 
along the horizontal axis comes from the fact that the calcula- 
tions are performed on finite-size clusters. The gray curves are 
the spin wave dispersion in the pure case. 

In order to consider the origin of the anomalies, we 
perform detailed calculations in one dimension. We take 
N — 256 and the random average for 500 configurations. 
Figure 3 shows the results for different values of x at 
Wimp/W = 0.5 {W = 2t in one dimension). The precise 
calculations in one dimension reveal the detailed struc- 
ture of the anomalies; an anti-crossing structure with a 
gap opening, shadow bands and a significant broaden- 
ing. These remarkable features become conspicuous for 
the values of x far from the quarter filling x = 0.5. 

An important point in Fig. 3 is that the anomalies ap- 
pear near the Fermi wave number fcp which is given by 
tt{1—x) in one dimension. This indicates that the anoma- 
lies come from a singularity of 2kp momentum transfer 
between q :^ kp and q ~ — fcp- There are also small 
anomalies at g ~ ir — kp and —Tr+kp. We note that these 
anomalies can be seen even in the case with a single im- 
purity. Namely, positions of anomalies in g-space scale 
with kp , while the intensities of the anomalous parts de- 
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Fig. 3. Spin excitation spectra in one dimension for (a) x = 0.5, 
(b) X = 0.4, (c) X = 0.3 and (d) x = 0.2. The arrow in each 
figure indicates the magnitude of kp. 

We speculate that randomness induces modulation of 
effective exchange coupling due to the Friedel oscillation. 
The random potential causes the 2fcF singularity in the 
charge channel, which appears in the effective ferromag- 
netic coupling through the magnon-electron interaction. 
Scattering of magnons by the oscillating components of 
the exchange couplings creates the anti-crossing struc- 
ture. The Friedel oscillation occurs universally in itiner- 
ant electron systems with randomness although its effect 
might be weaker in higher dimensions. Thus, the anoma- 
lies in three-dimensional systems in Fig. 2 may be an 
analog of the anti-crossing due to the Friedel oscillation. 

We discuss our findings in comparison with experi- 
mental results. Our results reproduce the anomalies of 
spin excitations reported in low-Tc compounds qualita- 
tively: (i) The spectrum shoY/ 's a sity nificarit broaden- 
ing near the zone boundaryMMEil'B (ii) The soften- 
ing near the zone boundaryHj'Ej'EB' can be interpreted 
as a consequence of the branching or the anti-crossing 
in our results. The lower-energy branch has a substan- 
tial weight and is easily identified in neutron scattering 
experiments, while the higher-energy one may be diffi- 
cult to identify due to a deficiency of energy resolution 
in esjperiments. (iii) The anti-crossing or the gap open- 
ingE2lE3' have been reproduced. Thus, the randomness 
is reasonably considered as a hidden parameter to give a 
qualitative description for the anomalous spin dynamics 
in low-Tc manganites. 

Several other scenarios have also been proposed to ex- 
plain the spin excitation spectrum in low-Tc mangan- 
ites. As an origin of the softening, some additional el- 
ements to the DE -mechanism have been examined sunh 
as magnon-orbitalES' and magnon-phonon, couplings.C^' 
Magnetic origin has also been proposedxB' For the ori- 
gin of the broadening, magnon-phononcJ' or magnon- 
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electron couplinga23) have been discussed. These argu- 
ments attribute the origin of the spectrum change by the 
A-site substitution to the modification of the bandwidth 
through some additional couplings to the spin degrees of 
freedom. As discussed previously, however, the actual 
changes of the bandwidth in these compounds are small. 
Therefore, it seems difficult for these scenarios to explain 
the large changes of the spin excitations in a quantitative 
manner. Moreover, neither the anti-crossing nor the gap 
opening has been reproduced by them. 

Our results give a comprehensive scenario to under- 
stand the systematic changes of spin excitations in the 
A-site substituted manganites, since the change of the 
randomness are expected to be large, as previously dis- 
cussed in reproducing the systematic decrease of Tq. We 
suggest that the randomness gives important effects on 
various thermodynamic properties of A-site substituted 
manganites. 

We now propose some experimental tests for our sce- 
nario. One is the sample dependence. Our scenario pre- 
dicts a correlation between the sample quality and the 
anomalies in the spin wave spectrum. Therefore, even in 
compounds with the same chemical formula, the extent 
of the anomalies would depend on the purity of sam- 
ples which is measured independently, for instance, by 
the residual resistivity. The second test is to investigate 
the location of the anomalies in the whole Brillouin zone, 
and to make a comparison with information of the Fermi 
surface, which is obtained by independent experiments, 
for instance, by the angle-resolved photoemission. The 
anomalies appear at the Fermi wave number in our sce- 
nario. Doping dependence is also crucial in this test. The 
third is to observe the high energy excitation near the 
band top q ~ (7r,7r,7r). In our results, there remains an 
excitation with substantial weight and relatively sharp 
peak-width near the band top even under strong ran- 
domness. This may be observed by the pulse-neutron 
scattering. All these experimental tests are feasible and 
may help to discuss the origin of the anomalies in the 
spin excitation as well as the role of the A-site substitu- 
tion in CMR manganites. 

Finally, we comment on higher-order corrections in the 
l/S expansion. The aim of this Letter is to study the 
systematic changes of the spectrum by the A-site substi- 
tution of manganites. In this case, the parameter \/S 
should be kept fixed, and the 1/5' expansion does not 
contribute to the spectrum changes. Furthermore, it has 
been shown that the contribution of the higher order 
terms are very small at x c^LL/lS, where we are interested 
in, at least in the pure caseEa-EB' Therefore, we conclude 
that the lowest order terms in the 1/5 expansion is suf- 
ficient to discuss the randomness effects as a first step, 
although quantitative estimates of the 1/5 corrections in 
the presence of randomness are left for further studies. 

To summarize, we have investigated effects of ran- 
domness on the spin excitation spectrum in the double- 
exchange model. The spin dynamics in the ground state 
is calculated by the spin wave approximation in the low- 
est order of the 1/5 expansion. Randomness induces 
some anomalies in the spectrum near the Fermi wave 
number such as the broadening, the anti-crossing and 



the gap opening. We speculate that the origin of the 
anomalies is modulation of effective exchange coupling 
by the Friedel oscillation. The anomalous spin dynamics 
in CMR manganites is reproduced qualitatively. Since 
the anomalies in experiments become conspicuous in A- 
site substituted compounds which have low Curie tem- 
perature, our results suggest that the randomness is an 
important element in the A-site substitution which has 
previously been understood as the bandwidth control. 
We have discussed our results in comparison with other 
theories and proposed several crucial experimental tests. 
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